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A RELATION BETWEEN BURNING VELOCITY AND QUENCHING DISTANCE 
By A. E. Potter, Jr. and A. L. Berlad 


SUMMARY 

A thermal quenching equation and the Semenov turning velocity equa- 
tion were combined to give ej^jressions for the product and quotient of 
burning velocity and quenching distance. 

The equation for the product was tested using published data. Experi- 
mental results for hydrocarbon-oxygen- inert gas, hydrocarbon -nitrous 
oxide, hydrogen-air, hydrogen-bromine, and the ethylene oxide deconposition 
flames were fairly well correlated by the product equation. Hydrocarbon- 
nitric oxide and ethylene oxide deconposltion data were not correlated. 

It was concluded that the product is proportional to a term which has the 
character of a thermal dlffusivity and a term which increases with the 
activation energy of the flame reaction. The quotient of btirning velocity 
and quenching distance is proportional to the average flame reaction rate 
and to terms mildly dependent on the flame reaction activation energy and 
the transport properties . Average and maximum reaction rates and maximum 
space heating rates were calculated for several flames xising the quotient 
equation and published data. The calculated maximum space heating rate 
for isooctane was about one and a half times the e3q>erimental value. 


INTRODUCTICW 


Burning velocity and quenching distance are closely related. Several 
hypotheses of flame quenching relate the two (refs. 1 to 3). For exanple, 
reference 1 gives the following equation derived from the thermal 
considerations : 


Ud 




(Tf - Ti)]V2 

^ T 


( 1 ) 


The product of burning velocity U am quenching distance d for a given 
tiibe geometry is directly proportional to the thermal dlffusivity a. and 
a term related to the chemical properties of the flame (the ignition tem- 
perature Ti is considered a chemical property). For flames with similar 
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thermal diffusivities and ignition temperatTires (hydrocarbon-air flames, 
e.g.) the product Ud should be a constant. Experiments show this to 
be true (ref. 4). However, if one wishes to use eguation (l) for com- 
parison of flames having different chemical properties, the problem of 
definition of the ignition tenperature arises. Equation (l) was derived 
from a slnple concept of flame quenching and the Mallard-IieChatelier 
btirning-veloclty equation. A better relation of TJ to d can be derived 
using more sophisticated quenching-distance and burning- velocity equations. 
Thiis, an equation for Ud was derived in reference. 3 by combining a 
dlffusional quenching-distance equation with a diffusional bvirning- 
veloclty equation. Only a limited test of the equation was made because 
of the lack of suitable data. In a similar way, the thermal quenching 
equation (ref. 5) can be combined with the Semenov biJirnlng-velocity equa- 
tion. The equations obtained in this way do not include an ignition 
tenperature. However, a tenperature Tq which is the minimum temperature 
to which a flame can be cooled before it is quenched does appear. This 
tenperature is not equivalent to the- usual experimental ignition tenpera- 
ture, which, is the tenperature at which unburned gas spontaneously 
ignites . 

Two combinations are interesting, the product Ud and the quotient 
U/d. In this paper, expressions for these combinations are discussed 
and tested using published experimental data. 


SYMB0I5 

A Arrhenitis pre -exponential term, (molecules/cc)^~“/sec 
a fuel concentration, molecules/cc 

b oxidant concentration, molecules/cc 

Cp specific heat at constant pressure, cal/(g)(deg) 

Cp average specific heat at constant pressure,^ cal/ (g) (deg) 
d quenching distance, cm 

E activation energy, cal/mole 

F fraction of heat prodToced by flame which must be retained in order 
for flame to exist 

f dimensionless geometrical constant 

G dimensionless geometrical constant 
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k thermal conductivity^ cal/(cm) (deg) 
ra apparent order of flame reaction 
n order of flame reaction with respect to oxidant 
R gas constant, cal/(mole) (deg) 

T temperat-ure , °K 

U burning velocity, cm/ sec 

w flame reaction rate, fuel molecules/(cc) (sec) 


w 


average flame reaction rate. 



fuel molecules/ (cc) (sec) 


a 


thermal diffusivity, k/pcp, cm^/sec 


(X* apparent thermal diffusivity. 


1/2 


^Po'^p,(l'=p,f/ 


n2/i 


cm‘="/ sec 


p density, g/cc 

Subscripts : 
f flame 

1 ignition 

Z lean limit 


o Initial 
q_ quenched 


THEORY 

Product of Quenching Distance and Burning Velocity 
The thermal qiienching equation of reference 5 can be written as 


PoCp,qWq^ 


( 2 ) 
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The subscript r, referring to the reaction zone, ^as previously used. 
Since the concept of reaction zone Is not strictly necessary for deriva.- 
tion of equation (2), the more general subscript q, referring to. the 
”about-to-be. quenched" flame, is used here. 

The Semenov burning-velocity equation (ref. 6) can be written as 




2kpW^ 

^'oPo^ 


( 3 ) 


Two terms which normally appear in this equation are not included here. 
These are the ratio of heat. .and mass diffusivitles ^and the ratio of moles 
of product to moles of reactant. These factors have been omitted from 
equation (3) because it is not known how they would appear in the quenching 
equation. This omission is quite jTistifiable for flames in which these 
factors are close to unity (this includes most flames), but must be kept 
in mind when discussing flames in which they differ greatly from unity. 


Equations (2) and (3) can be multiplied to give 


Ud = a* 



( 4 ) 


The term a"**” is called the apparent thermal diffiislvity and is defined 
in the symbol list. The term G is the geometrical constant. The 
factor F (fraction of the total heat of- combustion that must be retained 
by the flame for It-to exist) determines the temperature of the "aboutr 
to-be quenched" flame, which in turn determines the reaction rate 

in the "about -to-be quenched" flame. The value of Wg^ relative to* W|- 

is fixed not only by F, but also by the activation-energy of the reaction. 
For a given F value, a reaction having a small activation energy gives 
a smaller ratio w^/wg^ than one having a large activation energy. In 

general, the product Ud (for a given tube geometry) is directly prjopor- 
tlonal to a thermal diffuslvlty and to a term involving the chemical 
nature of- the flame. The close similarity of this equation to equation 
(l) is thus apparent. -i 

Calculation of apparent thermal dlffusivlty . -All the quantities 
which appear in the apparent thermal dlffusivlty are known or can be 

easily calculated except-the thermal conductivities". The assun5)tions 
and methods used for their calculation are given in' the appendix. For 
consistency with the original derivation of the th ermal quenching equa- 
tion (ref. 5), the heit capacity and thermal conductivity in the quenching 
eqaiation both refer to unburned gas at 0.7 flame tengjeratur e . Becavise 
of the vague plct\3re used for that derivat-ion, it is by no means certain 
that these are referred to the correct conditions . Only a better 
understanding of the flame„quenching can decide this question. 
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Calculation of Fw^/wg^. - In order to calculate Fw^/w^, it is 

assumed that Semenov's formulation of the average reaction rate is adequate 
for con^iuting w, and that it applies to both the free and quenched flame. 
Semenov's formulation of the average reaction rate for stoichiometric 
flames for zero, first, and second order reactions can he obtained from 
reference 6. Generalizing these to the case of a fuel-oxygen reaction 
of m"^^ order gives 


^f = 




n 


RToTf 


im 


RT‘ 


E 

RTf 


E(Tf - To) 


E(Tf - Tq) 


( 5 ) 


This generalization was made in order to permit the use of a fractional 
reaction order if necessary. Eq.uation (5) is correct for zero, first, 
and second orders, and consequently can be expected to be approximately 
correct for fractional orders close to zero, one, and two. 

From reference 5, (assuming that is the average flame temperature 

in the "about -to -be quenched" flame, and that heat capacity is constant) 

P = (T<^ - To)/(Tf - Tq) (6) 

Combining equations (5) and (6) and an expression for Wg^ analogous to 
equation (5), and assuming A to be Independent of tenperature. 



From equation (7), it can be seen that values of the flame activa- 
tion energy are required to calculate Pw^/wg^. Now the activation energy 

obtained from biirnlng-veloclty data depends somewhat on the assumptions 
made concerning the reaction order and the temperature dependence of the 
Arrhenius constant. The best procedure appears to be to use a consistent 
set of activation energies. Assuming the Arrhenius constant to be tem- 
perature independent (as done in the derivation of eq. (?)) and a bl- 
molecular reaction, the author of reference 7 fo\md that the low- 
temperatm'e activation energy for propane oxidation of 38 kilocalories 
per mole was a satisfactory value for propane-air flames. The authors of 
reference 8 found that the activation energy was nearly the same for most 
hydrocarbon-air flames. On the basis of these results, an activation 
energy of 38 kilocalories per mole was used for all hydro carbon- oxygen- 
inert gas flames. For some of the other flames, no activation energies 
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are available. However, the activation energy for hydrocarbon-air flames 
is directly proportional to the lean-limit flame ten^jeratijre (ref. 8). 
Activation energies were calculated by this method for hydrog-en-air , 
hydrocarbon-nitric oxide, and hydrocarbon-nitrous oxide flames, using 
lean-limit flame temperatiu'es for hydrogen-air, pentane -nitric oxide, and 
methane -nitrous oxide flames. 'The lean-limit flame ten^peratures were 
multiplied by 24 to make them consistent with a v^ue of 38 kilocalories 
per mole for propane-air flames. Low-temperature activation energies 
were \:iBed for the hydrogen -bromine and ethylene oxide decomposition 
flames. 

Another quantity necessary to calculate Pwf/wq is the reaction 
order m. Values of m are not available for most flames. Fortunately, 
Fw^/wq^ is quite Insensiti've^to the value chosen for m. Since the order 

is close to two for most of- the flames discussed herein, the order was 
taken to be two for all flames except the ethylene oxide decomposition 
flame, for which an order of one was assumed. 

The minimum temperature to. which a flame can be cooled and continue 
to exist is Tq. In order to determine Tq from equation (6), the value 
of F must be known. For lack of any other information, it will be 
assimied that F is the same for all flames. It is found that this 
assumption- is a good one. Using the set of activation energies just dis- 
cussed, a suitable value (found by trial and error comparison of experiment' 
and theory) is 0.78. If the activation energies afe decreased by 30 
percent to bring them In line with the values given in reference 8, the 
best value for F drops to about 0.70. 

Some values of (Fw^/w^)^/^ for several different-flames along with 
the activation energies used are given in table I. 

Once a value for F has been chosen (0.78 in this case), the size 
of (Fwp/Wq)V2 is determined solely by the size- of the activation energy 

relative to the flame temperattire . This is clearly shown in table I. 

For flames with similar activation energies and flame temperatures, very 
little difference is seen in the (F?f^/?7q)V2 values. For these flames, 

one can -write equation (4) as 


Ud oc a* (8) 

Sgurces -of experimental data . - The values of quenching distance and 
burning velocity used for a test of equation (4) are shown in table II. 
Quenching distance is fairly reproducible, so that quenching data were 
taken directly from piiblished results. All quenching distances are for 
plane parallel plates. When necessary, quenching data for cylinders 
were converted to plane parallel pla-fces by multiplying them by 0.66 





IT* 


NACA TN 3882 


7 


(ref. 9). Burning velocities vary according to the method of measxirement , 
so the data were arbitrarily adjtisted to a common basis whenever possible. 
Sets of data were chosen which contained a value for propane -air flames. 
These data were multiplied by a factor to make the stoichiometric propane 
burning velocity equal to 40 centimeters per second at 25° C. Data from 
references 10 and 11 were adjxisted to a common basis through the burning 
velocities given for benzene-air flames. Stoichiometric burning veloc- 
ities were used when available; otherwise, mayiTmim values were used. 
Tenperature exponents of burning velocity were assumed to be 1.7 for 
hexane and hexene-1 flames, and to be 2 for hydrogen and propane with 
argon and helium "air". For hydrogen-bromine flames, the btirnlng velocity 
on the 1,5-centimeter burner was used rather than that for the 1.0- 
centlmeter burner. The values used for the quenching distance and burning 
velocity of ethylene oxide are believed to be the best available. 

Comparison of theory and experiment, - The test of equation (4) is 
shown in figure 1, where the observed product Ud is plotted against 
the calculated product for the 30 flames given in table H. Generally, 
hydrocarbon-oxygen-lnert gas flames are well correlated, although flames 
at elevated initial temperatures show a trend away from a mean line. It 
seems likely that the assumption of a constant F value is incorrect for 
these flames. The significance of this is not clear at this time. 

Several other flame types are well correlated; these are hydrocarbon- 
nitroios oxide, hydrogen-air, and hydrogen-bromine flames. The hydrocarbon- 
nitric oxide and ethylene oxide decomposition flames are not correlated. 
These flames present considerable experimental difficulty, and this may 
explain the discrepancy. 

From these results, one can conclude that the product Ud is directly 
proportional to a term having the nature of a thermal dlffiislvity and a 
term that increases with increasing activation energy of the flame 
reaction. 

'pypirlcal relation between U and d . - Occasionally, the relation 
between U and d for chemically similar flames, such as hydrocarbon- 
oxygen-nitrogen flames, is of interest. According to equation (8), for 
systems such as this, the product is directly proportional to the apparent 
thermal dlffusivlty. From the definition of the apparent thermal dlffusi- 
vity, it can be seen that this quantity varies directly as the initial tem- 
perature, inversely as the pressure, and directly as flame temperature to 
a power describing the temperature dependence of the thermal conductivities 
and heat capacities. By comparison with experiment, it is found that a good 
value for the flame temperature exponent is 1. A plot of Ud for 
hydrocarbon-oxygen-nitrogen flames at 1 atmosphere against T^T^ is shown 

in figure 2. By least squares, it is found that the relation has the 
following form; 


Ud = 1.204xlO"^ToTf 


(9) 
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Quotient of Btiriilng Velocity and Quenching Distance 
By dividing equation (3) hy equation (4) and rearranging, one obtains 



The terms under the square root sign do not vary much from flame to 
flame. The reciprocal of— the first of these, Fw^/wg^, is discussed 

previously, and it was found to be a mild function; of activation energy. 
The other term contains the ratios of thermal conductivity and heat 
capacity in the free and "about-to-be quenched” flame. This term can be 
expected to be nearly constant. Its value depends largely on the assunp- 
tlons made about the tenperature and composition of the "about-to-be 
quenched" flame. Thus, for similar flames (e.g., hydrocarbon-oxygen- 
inert gas flames). 


U '^f 

^ oc 


Using the assumptions previously made concerning Ftr^/Wq^, It Is 

possible to use equation (lo) to calculate average reaction rates from 
experimental burning velocities and quenching distances. This has been 
done for some representative flames, and the results are shown in table 
III. The rates are quite sensitive to the choice of — F; a 3-percent 
change in F changes the calculated rate about 10 percent. To compare 
the flames on a common basis, the average rate calculated from the data 
was corrected to an initial tenperature of 25° C -using equation (5). 

Two values of particular interest are those for pr^ane with argon "air" 
and with helium "air". Since replacement of argon by hellttm does not 
affect flame temperatvire , it Is often ass-umed that reaction rates in the 
two flames are identical. The reaction rates calculated from equation 
(lO) differ considerably for the two flames. Since two experimental 
results are combined in the quotient, the dlffei:ence may be experlmentaL 
error. Other evidence that a real difference may exist is the fact that 
the pressiire exponents of quenching distance (related to the order of 
the reaction) are different for the two flames (ref. 12). This seeming 
difference could mean that an Important step in the- complex chain reaction 
of the coaibustion process is a three-body reaction involving the inert 
gas. . — 



The maxim-urn reaction rate attained in the flame is frequently of 
interest. Provided the flame reaction can be approximated by a simple 
one-step reaction of the hind treated by Semenov, the maximum rate can 
be found from -the average rate. Again generalizing from equations in 
reference 6, the reaction rate in the flame is 
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( 12 ) 


The temperattire corresponding to the maximum rate can he obtained by- 
differentiation of equation (12) with respect to ter^ieiature and setting 
the result equal to zero. This optimum temperature is 


"^opt “ 


^ 

mRTj, 


(13) 


To find the relation of the ms,ximum reaction rate to -the average rate, 
equations (5), (12), and (13) are confeined. Then, 


w. 


'max 


— E (^f ~ ^^q) 

— TJ’^ ■ ■■■ 


= Wf. 

^ E 


upi 0-m 


(14) 


Values for for some representative flames are given in table m, 

along -with maximum space heating rates (derived from v-ma-r multiplying 

it by the heat of combiistion of the ftiel). These maximum space heating 
rates maybe valuable for cor^jarison of different fuels. They are not 
absolute values because of the many approximations Introduced. By way 
of co3i5)arison, LongwelL’s val-ue of -fahe space heating rate of isooctane- 
air flames at an initial ten^jerat-ure of 300° K and 1 atmosphere is 
2.1KL0® B-tu per cubic foot per ho-ur (ref. 13), and the calculated value 
for the same initial conditions is 3.6x10® Btu per cubic foot per hour. 

The values of reaction order m used are given in table HI. They were 
calculated from -the pressure dependence of quenching distance using the 
procedures of reference 14, and are corrected for the effect of pressure 
on flame ter^ierature , The references cited in table HI are the so-urces 
of the pressure exponents. 


CCWCLUSIOfRS ^ 

O 

Combination of a thermal quenching equation 

1/2 


and a b\3rning velocity equation 

B-uming velocity oc [(Reaction rate) (Transport property)] 
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gives the following expressions for the product and quotient of burning 
velocity and quenching distance: The product is 

(Burning velocity) (Quenching distance) oc Transport property 



In practice these statements must be modified somewhat, since reaction 
rates and transport properties differ in the free and quenched flames. 
However, if similar flames (e.g. , hydrocarbon-oxygen-lnert gas flames) 
are compared, the proportionalities are correct as stated. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 15, 1956 
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APPENDIX - CALCULATION OP FLAME TEMPERATUEE, THERMAL 
CONDUCTIVITY AND HEAT CAPACITY 

Flame tesmperatiires for tlie hydrocarLon-alr flames at an intlal tem- 
perature of 378° K were calculated from the tables given in reference 15. 
All other flame teii 5 )eratures are values taken from the references in 
which burning velocity is given. 

Thermal conductivities for the Lennard- Jones (6-12) potential were 
calculated for individual components of the mixtures for which conduc- 
tivity values were desired, using the procedures and data given in 
chapter 8 of reference 16. Mixture conductivities were calculated by 
the methods described in reference 17. The viscosity equation given by 
Hirschf elder (ref. 16) was used in the Llndsay-Bromley mixing rule rather 
than the Sutherland equation. For mixtures containing water vapor, the 
procedure given on page 600 of reference 16 was used to , obtain force con- 
stants for use in Hirschf elder ' s viscosity equation. 

For kq_ values, the conductivity of the unbumed mixture at 0.7 Tf 
was calculated. Since heat is lost to the wall through unbumed mixture, 
this seems q.ualitatlvely reasonable. 

For kf values, the conductivity of the blamed gases at T^ was 
calculated. For simplicity, dissociation was ignored, and k^ was cal- 
culated for the simple products of reaction at Tp. 

Heat capacities were calcoilated from, the tables of reference 18. 

For the heat capacity of the unbumed mixture at 0.7 T^ was cal- 

culated. For Cp^f, the average heat capacity of the undissociated 
products from Initial temperature to flame temperature was calculated. 
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TABLE I. - VALUES OP (EWf/wg^)V2 FOR SEVERAL FLAMES 


Flame system 

(Fw^/W(^)V2 

Activation 

energy, 

kcal/mole 

Flame 

ten^jera- 

ture, 

Tf, °K 

Propane -alr^ 

2.33 

38 

2289 

Isooctane -air 

2.32 



2302 

Methane -oxygen 

2.06 



3020 

Hydrogen-air 

1.69 

23 

2380 

Ethane -nitrous oxide 

1.80 

29 

2740 

Ethane -nitric oxide 

2.69 

49 

2840 

Hydrogen-hromlne^ 

4.32 

40.2 

1490 

Ethylene oxide decomposition 

8.56 

52.7 

1217 


^•Initial temperature Tq, 378° K, 
^Equivalence ratio, 2/3. 




H 

TABii; II. - HtnaoBB VELocinr, cpncuHa diSeabck, a® AcrwAnos hsikjt ^ 

FOR VARI003 STOtCHIOIffiroiC FIAMB8 


Flame aystem 

Enitlal 
teag>er- 
ature, 
Po, °K 

Burning 
velocity, 
U, cm/aec 

Refer- 

ence 

for 

burning 

velocity 

Quenching 
distance, 
d, cm 

Reference 

for 

gpenchlng 

distance 

Activa- 
tion 
energy, 
E, heal 

Refer- 

ence 

for 

activa- 

tion 

energy 

Flame 

teoQ>era- 

ture, 

■Ef, ^ 

Prcqiane-alr 

S78 

®58.3 



0.184 

r 1 

4 

38 



2289 

n-Hexane-alr 



“57.3 

10 

.182 



^38 



2300 

n-De cane -air 



»55.9 

11 

.183 







2303 

iBooctane-aJx 



»44.5 

19 

.233 







2302 

Hexene -1-air 



^.1 

10 

.166 







2324 

Decene-l-alr 



“60. 6 

11 

.175 







2322 

Bensane-alr 



“63.6 



.171 







2365 

n-Butyl bensene-alr 



“53.5 



.203 







2559 

Propane-17S{ oxygen-nitrogen 

313 

“20.3 

1? 

.377 

20 





2051 

Erqpana-2lJ{ oxygen-nitrogen 



“42.3 



.199 







2253 

Prqpane-SO^ oxygen-nitrogen 



“91.8 



.087 







2558 

Parcrpane-SO^ oxygen-nitrogen 



“202 



.051 







2844 

Prqpatie-70)t oxygen-nitrogen 



“312 



.032 







2970 

Protpane-argon "air" 



“90.1 

21 

.104. 

12 





2557 

Propana-hellum "air" 



“158 

' 


.253 

12 


' 



2557 

PrcipanB-alr 

400 

“64.1 

7 

.183 

22 


f 



2313 

Propane-air 

483 

“69.0 



.159 




i 



2351 

Propane-air 

5B8 

“115.5 



.144 




j 



2384 

Ethylene -air 

298 

“69.9 

10 

.125 


3 





2362 

Methane-air 



“32.8 

23 

.25 

24 





2214 

Msthsne-axygen 



“344 ■ 



! .03 







1 3020 

Hyflrogen-alr 

313 

“216 



.0495 

14 

“24 

8 

2380 

Acetylene-nitrous oxide 

298 

160 



.056 


4 

C29, 

25 

2940 

Ethyiene-nltrons oxide 



120 



.048 




. 



2820 

Ethane-nitrous oxide 



90 



.059 





' 


2740 

Acetylene -nitric aside 



85 



.24 



“4 

9 

26 

3080 

Ethylene -nitric oxide 



40 



.44 







2920 

Ethane -nitric aside 



23 



.86 







2840 

Hydrogen (60^) -broolue (40iO 

323 

28 

2- 

1 

.59 

27 

40.2 

20 

1490 

Ethylene oxide decoigiosltioD 

298 

2.7 

29 

4.6 

30 

52 

.7 

29 

1217 


®7aliiea edjusted aa described In the text. 

^Estimated. 

; *^Calculated Aron, lean- Unit flaiae teaperatnre B “ 24Tj 
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TABLE III. - REACTION RATES AND MAXIMLIM SPACE HEATING RATES 


Flame system 

Initial 
temper - 
at'ure. 
To, °C 

Order of 
reaction, 
m 

Average 
reaction 
rate, w, 
moles/sec 

Maximum 

reaction 

rate, 

Wmax, 

moles/sec 

Maximum space 
heating rate at 
1 atm, 

Btu/ (cu ft) (hr) 

Hydrogen -air 

25 

s-2.17 

169. 

463, 

108X100 

Methane -air 



^1.60 

1.60 

5.07 

3.9X100 

Ethylene -air 



c^l.60 

5.00 

15.1 

19.3X100 

Propane -air 



ai.60 

1.04 

3.27 

6.44X100 

Propane - argon 







"air" 



^1.71 

3.82 

11.3 

22.4X100 

Propane - hellnm 







"air" 



^1.63 

2.94 

8.40 

16.6X100 

Isooctane-air 



ai.52 

.24 

.72 

3.6X100 


s-Ref. 14. 
^Estimated. 
^Ref. 3. 
<iRef. 12. 
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EACA TN 3882 
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Figure 1. - Comparison of experimental and calculated values 
of 'burning-velocity - quenching -distance product Ud. 
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Figure 2. - Relation of empirical buming-velocity - q,uenching-dlBtance product 
Ud to empirical temperature product for hydrocarton-oxygen-nitrogen 

flames . 


NACA - Langley Field, Va. 



